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Starting from a uniform d-wave superconducting phase we study the 
energy cost due to imposed unidirectional defects with a vanishing pair- 
ing amplitude. Both renormalized mean-field theory and variational Monte 
Carlo calculations within the t-J model yield that the energies of inhomo- 
geneous and uniform phases are very close to each other. This suggests 
that small perturbations in the microscopic Hamiltonian, might lead to in- 
homogeneous superconducting phases in real materials as observed in recent 
scanning turmeling microscopy on Ca2-iNaa:Cu02Cl2. 

PACS numbers: 74.72.-h, 74.20.Mn, 74.81.-g, 75.40.Mg 

1. Introduction 

Recent progress in spectroscopic techniques has provided a wide variety 
of interesting data concerning electronic states of the high- Tc superconductors. 
For example, scanning tunneling microscopy (STM) on different cuprate fami- 
lies Ca2-a;Naa;Cu02Cl2 and Bi2Sr2Dyo.2Cao.8Cu208+5, has revealed short-range 
unidirectional charge domains coexisting with inhomogeneous c?-wave supercon- 
ductivity T^. In particular, it has been found, that the doped holes primarily 
enter oxygen sites leading to a bond- centered charge pattern with a period of four 
lattice spacings. Motivated by this result, we have recently shown [2J, that such 
a charge order might be naturally interpreted in terms of a valence bond crystal 
[3], i.e., paramagnetic phase with both spatially varying bond charge hopping 
and short-range antiferromagnetic (AF) correlations. In this case, an inhomoge- 
neous antiphase domain resonating valence bond (ttDRVB) phase was obtained 
by assuming a 7r-phase shift in the superconducting (SC) order parameter across 
domain walls (DWs). While the antiphase solution is particularly intriguing since, 
in contrast to its inphase counterpart, offers a simple explanation of the suppres- 
sion of the effective interlayer Josephson coupling observed in some stripe-ordered 
high-Tc compounds [HIS], both types of the modulation of the SC order parameter 
are the subject of intense ongoing studies [3 [71 [HI ISj ■ Therefore, in this paper 
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we shall study the energy cost due to imposed defects with a vanishing pairing 
amplitude (no a 7r-shift is assumed across the DWs) and compare the resulting 
charge modulation with the corresponding one found in the ttDRVB phase. 

2. Model and the approach 

We investigate a t-J model Hamiltonian, 

^ = E (4s. + h.c.) + jE • (1) 

where = (1 — n,;^_(j)c|^ is the Gutzwiller projected electron operator and 
use a renormalized mean field theory (RMFT) in which the local constraints of 
no doubly occupied sites are replaced by statistical Gutzwiller weights gl^ idij) 
for hopping (superexchange) processes, respectively [TD]. Hence the mean-field 
Hamiltonian reads, 

Hmf = E 5li(c,>S\<T + ^^-^O - ^^^n,^<y 

with the Bogoliubov-de Gennes self-consistency conditions for the bond- Xji = 
i'^],a^i,cr) pair-order A^^ — {Cj_„c^^^) = (c^ _^Cj^) parameters in the unpro- 
jected state. We consider here the so-called modified Gutzwiller factors, 

J ^ A{1 ~ nhi){l - Uhj) 

+ 8nH^nhJp-J (2) + 16/3+ (4) ' 



a,,+%{l-nh^nhi)\x^JV + l&\x^Jr ^' 

where = {'^ - nl-}{l ~ nlj), Pf^iji) = |Ay |" ± \xij\^ while are local hole 
densities. By including the effects of the nearest-neighbor correlations Xij and 
Aij they are known to give a better agreement with a more accurate Variational 
Monte Carlo (VMC) technique [2]. Hereafter, we shall assume a typical value 
t/J = 3 and fix the doping level x = 1/8. Finally, using unit cell translation 
symmetry [11], RMFT calculations were carried out on large 256 x 256 clusters 
at a low temperature 13 J = 500 approaching thermodynamic limit. 

3. Results and discussion 

In Fig. [T] we show the hole profiles as well as the values of the bond- and 
pair-order parameters across the unit cell found in the ttDRVB (top) and inphase 
DRVB (bottom) state. The obtained modulations clearly reflect the competi- 
tion between the superexchange energy Ej and kinetic energy Et of doped holes. 
However, a detailed charge profile depends on the assumed type of the SC order 
parameter. On the one hand, suppression of the pair-order amplitude Ay along 
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Fig. 1. (a,b) Hole density n^i and variational parameters: (c,d) Ai^i+a as well as (e,f) 
Xi,i+a found in the ttDRVB (top) and DRVB (bottom) phase. Solid (open) circles in 
panels (c-f) correspond to the x (y) direction, respectively. 

the DWs automatically involves a deviation of the bond-order parameter Xij from 
the value found in the areas with finite Ay . Remarkably, the deviation is par- 
ticularly strong in the case of the antiphase SC order parameter. On the other 
hand, the absence of the tt shift across the stripe boundary in the DRVB phase 
allows the system (as confirmed by the VMC method [8J) to avoid a reduction of 
Aij on the adjacent vertical bonds which remains almost intact. Therefore, the 
charge redistributes from the hole rich areas with enhanced Ay in the ttDRVB 
phase towards DWs with vanishing A^ in the DRVB state (see Fig.[T]). 

In order to appreciate better the reason of a different charge profile in both 
phases we show in Fig. [2lja-d) the corresponding short-range AF correlations, 

ST = ^\9i^+J\X^,^+.? + \^^,^+c.?), (5) 

with a = {x^y}, as well as bond charge hopping, 

^" = 2g[,+„i?e{xM+a}, (6) 

across the unit cell. Here one finds that a local reduction of the SC order param- 
eter (and the concomitant strong suppression of the superexchange energy on the 
related bonds) enables, in the ttDRVB phase, a large bond charge hopping along 
the DWs as in the usual stripe scenario [2]. In fact, it also determines the actual 
hole profile arranged in the way which minimizes the loss of the superexchange 
energy at the DWs. This can be easily accomplished by expelling the holes and 
strengthening locally the corresponding gl^ factors. In contrast, small modulation 
of Xij ill the DRVB phase results in a much weaker, with respect to the ttDRVB 
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Fig. 2. (a,b) Spin correlation 5f , (c,d) bond charge 7^", and (e,f) SC order parameter 
Af^ found in the ttDRVB (top) and DRVB (bottom) phase. Solid (open) circles corre- 
spond to the X (y) direction, respectively; solid line in panels (e,f) depicts the SC order 
parameter in the uniform d-wave RVB phase. 

one, modulation of both the spin correlations and bond-charge hopping. Conse- 
quently, the system does not have to further improve the superexchange energy 
at the defect lines but it rather tries to regain some kinetic energy released on the 
broken RVB bonds. This is reached by adjusting the hole profile and attracting 
the holes to the DWs which enlarges locally renormalization factors gjj . As a 
result, the DRVB phase has a very good kinetic energy being even slightly better 
than that of the uniform d-wave RVB phase (see Table I). Let us point out, 
however, that even though both the RMFT and VMC methods predict exactly 
the same hole profiles in the DRVB phase (as well as its remarkably good energy), 
a discrepancy appears concerning kinetic energy gain at the DWs, strongly en- 
hanced in the VMC method [5]. The difference simply follows from the fact that 
in the RMFT both the short-range AF correlations and bond-charge hopping are 
(X Xij ■ Hence its suppression involves a reduction of both the energy contributions 
unless the system is disposed towards a strong phase separation so that they can 
be further modified by the Gutzwillcr factors [T2] . 

TABLE I 

RMFT kinetic energy Et, magnetic energy Ej, and free energy F as well as VMC energy 
Eymc of the locally stable phases: ttDRVB, DRVB, and d-wave RVB one at a; = 1/8. 



phase 


Et/J 


Ej/J 


F/J 


Eymc/ J 


ttDRVB 


-0.8719 


-0.4518 


-1.3237 


-1.3359 


DRVB 


-0.8871 


-0.4662 


-1.3533 


-1.3647 


RVB 


-0.8863 


-0.4784 


-1.3647 


-1.3669 
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Finally in order to discuss the SC properties of our inhomogeneous phases 
we plot in Fig. [2l^e,f ) the modulus of SC order parameter, 

AL'' =5m+„|A,.,+„|, (7) 

across the unit cell. One of the key qualitative differences between the ttDRVB 
and its inphase counterpart is evident in this figure. Namely, while the SC order 
parameter deviates, in the regions between defect lines, only slightly in both 
states from the value found in the uniform d-wave RVB phase, the absence of the 
TT shift across the stripe boundary in the DRVB phase decouples the horizontal 
and vertical bonds constituting DWs. Therefore, in contrast to the ttDRVB phase, 
the latter retain the value of the SC order parameter of the uniform state. 

4. Summary and conclusions 

In this paper we have studied two possible modulations of the SC order 
parameter across the DWs: inphase and antiphase. Remarkably, we have found 
that the energy of the unidirectional modulated phases (especially of the inphase 
configuration) approaches the energy of the uniform d-wave RVB superconductor. 
In fact, the energy difference might be further reduced by the tetragonal lattice 
distortion that often appears in the high- Tc compounds [8] . We conclude therefore 
that the d-wave RVB phase is capable of efficient minimizing the energy cost due 
to unidirectional defects with broken RVB bonds which in turn might induce the 
charge modulation similar to that observed in the STM experiments [1] . 
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